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ABSTRACT 

The reionization of cosmic hydrogen, left over from the big bang, increased its 
temperature to > 10 4 K. This photo-heating resulted in an increase of the minimum 
mass of galaxies and hence a suppression of the cosmic star formation rate. The affected 
population of dwarf galaxies included the progenitors of massive galaxies that formed 
later. We show that a massive galaxy at a redshift z > 6 should show a double-peaked 
star formation history marked by a clear break. This break reflects the suppression 
signature from reionization of the region in which the galaxy was assembled. Since 
massive galaxies originate in overdense regions where cosmic evolution is accelerated, 
their environment reionizes earlier than the rest of the universe. For a galaxy of ~ 
10 12 M© in stars at a redshift of z ~ 6.5, the star formation rate should typically 
be suppressed at a redshift z > 10 since the rest of the universe is known to have 
reionized by z > 6.5. Indeed, this is inferred to be the case for HUDF-JD2, a massive 
galaxy which is potentially at z ~ 6.5 but is inferred to have formed the bulk of its 
3 x lO n M in stars at z > 9. 
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1 INTRODUCTION 

Following the epoch of cosmic recombination at z ~ 10 3 , 
neutral hydrogen pervaded the universe unti l it wa s reion- 
ized by the first galaxies dBarkana fc Loebl l200lf) . There 
has been much interest recently in developing observa- 
tional techniques that will reveal the redshift of reioniza- 
tion, z rc i. Existing methods provide conflicting constraints. 
The large-scale polarization anisotropies of the cosmic mi- 
crowave ba ckground (CMB) i mply 11 < z r < 30 at 95% 
confidence (|Kogut et all 12003! ) , whil e the size of the ion- 
ized regions around quasar s at z ~ 6 dWvithe fc Loebll2004 
iMesinger fc Haimanl 12004 IWvithe. Loeb, fc Carillil 120051 ^ 
imply a substantial neutral hydrogen fraction at that late 
time. Here we propose a new probe of z r , namely the im- 
print of reionization on the star formation histories (SFHs) 
of individual massive galaxies at z > 6. 

The photo-ionization of hydrogen raised its tempera- 
ture to > 10 4 K. This heating process resulted in an in- 
crease of the minimum mass of galaxies and hence a sup- 
pression of the cosmic star formation rate (SFR). We pre- 



* E-mail: 
(AL) 



barkana@wise.tau.ac.il (RB); aloeb@cfa.harvard.edu 



viously suggested (iBarkana fc Loebl l2000h that this global 
signature could be observed by detecting the corresponding 
fall in the number counts of faint galaxies. However, we pre- 
dicted that even the future James Webb Space Telescope 
could individually detect only a fraction of the affected low- 
mass galaxies. While the same signature may be more eas- 
ily detectable through the evolution of the supernova rate 
l|Mesinger. Johnson fc Ha iman 2005;), there is a generic rea- 
son that makes it advantageous to search for localized rather 
than global signatures of reionization. At high redshift, large 
fluctuations in the number density of galaxies in various re- 
gions causes reionization to occur at different times in differ- 
ent regions (|Barkana fc Loeb! 2004). Thus, a region of size 
< 10 comoving Mpc should show evidence for a relatively 
rapid reionization, while any global signature will be further 
smoothed out over a redshift interval of ~ 2-4. 

Since a 10 12 Mq halo forms out of matter in a region of 
a few Mpc in size, the SFH of an individual massive galaxy 
at high redshift should show a break in its SFR. This break 
reflects the suppression signature from reionization of the 
region from where the galaxy was assembled. We show that 
even a rather massive galaxy at z ~ 6 should show this 
break, since most of its mass is expected to have been in low- 
mass progenitors when the surrounding region was reionized. 
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Moreover, if this signature can be detected in multiple galax- 
ies which lie in regions that reionized at different times, then 
the spatial variation of reionization can be probed with this 
method. 

Massive galaxies are routinely being detected at red- 
shifts z > 6 jBouwens et~aH 120051 ; IStark fc Ellisll2006h . in 
some cases with sufficient d ata for a detailed analy sis of the 
stellar content. For example, Mobashcr et all 1120051 ) have re- 
cently reported the discovery of HUDF-JD2, a galaxy which 
is potentially at z ~ 6.5, although spectral lines were not 
detected and thus the redshift is based on the shape of the 
spectrum. If the high redshift is a correct inference, then the 
galaxy possesses almost 1O 12 M0 in stars. Further analysis 
of the spectral shape, in particular the Balmer break signa- 
ture of an established stellar population, indicates that the 
galaxy formed the bulk of its stars at z > 9 and could have 
reion ized it s surrounding regi on at z > 10 jPanagia et al.l 
2005). Also. lEvles et all (|2005l) analyzed two galaxies at z ~ 
5.8 th at had been found by Stanwav. Bunker, fc McMahonl 
(2003). They showed that these galaxies had stellar masses 
of > IO 1O M0 and had formed most of their stars at z > 7.5 
whil e a minori t y form ed closer to the detected redshift [see 
also lYan et ail i|2005l )]. 

It is therefore particularly timely to examine the ques- 
tion of whether the above-mentioned data imply that reion- 
ization occurred early in the universe. In this paper we con- 
sider the imprint of reionization on the SFHs of individ- 
ual massive galaxies like HUDF-JD2. Since massive galaxies 
originate in overdense regions where cosmic evolution is ac- 
celerated, we show that their environment reionizes earlier 
than the rest of the universe. In § 2 we present our for- 
malism for calculating the SFH within galaxies and show 
general results in § 3. We apply these results to HUDF-JD2 
and the other galaxies in § 4 and give our main conclusions. 
Throughout the paper, we assume a flat universe with cos- 
mological parameters h — 0.72 and fi m = 0.27 (which in- 
cludes an fi;, = 0.046), w ith the remaining en ergy density in 
a cosmological constant l|Spergel et al.ll2003l ). 



2 MODELING STAR FORMATION 
HISTORIES 

In order to model the SFH of a typical galaxy, we must ac- 
count for a number of effects that determine the availability 
of gas for star formation. We first consider the ability of gas 
in general to cool, and the ability of pressure to resist grav- 
ity - particularly after the gas is heated during reionization. 
Having analyzed these effects on the gas, we begin for each 
galaxy with the average accretion history of the dark matter 
halo that will eventually host it, and then calculate the ad- 
ditional effects that modify the accretion history of the gas 
relative to the dark matter. We then consider how the reion- 
ization history differs from region to region, and how this is 
correlated with the SFHs of galaxies in different regions. 

We begin by considering the gas in galactic halos. Star 
formation requires the gas to cool and collapse to high den- 
sity within a halo. We assume that galaxies form in halos 
above some minimum mass M coo i that depends on the effi- 
ciency of atomic and molecular cooling. A minimum virial 
temperature of ~ 1 4 K is required for significant ionization 
and thus efficient atomic cooling, and this imposes a mini- 



mum halo circular velocity of 16.5 km/s. On the other hand, 
if molecular hydrogen is present, then T v j r ~ 1 3 K is needed 
to excite its rotational and vibrational transitions, which 
yields V c = 3.7 km/s for molecular hydrogen cooling. Since 
feedback from stellar radiation m ay efficiently dissociate 
the H 2 in the intergalactic medium l|Haiman. Rees. fc Loebl 
Il997t ). we consider two possible cases, where either molec- 
ular cooling is effective or is not (atomic cooling sets the 
threshold in the latter case). In either case, once the min- 
imum V c has been determined, the minimum halo mass in 
which gas can cool efficiently is then 



M, 



0.14 



•(1) 



In addition, reionization heats the IGM, increasing pres- 
sure gradients dramatically and preven t ing gas from ac- 
cumulating in small ha los 



and p reventing ga s tro 
(iReesI 1 19861; lEfstathiou 



1992 



1996 



Thoul fc Weinberg! [f996l; iQuinn. Katz. fc Efstathiot 
Weinberg, Hernquist, & Katz 1 19971 ; |N avarro fc Steinmetzl 
19971 ; iKitavama fc Ikeuchil boOoT iDiikstra et all 120041 ). 

Specifically, we assume that after reionization gas contin- 
ues to accrete only into halos above a characteristic infall 
mass. For a small density fluctuation relative to the cosmic 
mean density, gravity overcomes pressure in a region above 
the Jeans mass 
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where T is the temperature of gas at the cosmic density, 7 
(which lies after reionization in the range 1-1.5) is its adi- 
abatic index, and {i is the mean molecular weight in units 
of the proton mass. However, a halo in the process of for- 
mation is significantly denser than the cosmic mean, so that 
even halos well below the Jeans mass can accrete some gas, 
especially at high redshift. Based on the above-mentioned 
numerical simulations, we adopt a minimum infall mass of 
Mj/4, which corresponds to a circular velocity of ~ 35 km/s. 
Note that it would be incorrect to use the filtering mass 
IIGnedin fc Hu¥l998) as the cutoff mass, since the filtering 
mass describes the cumulative time-averaged effect of pres- 
sure suppression while our results depend on the instanta- 
neous rate of gas infall into galactic halos. 

T he Jeans mass depends on the thermal ev olution of the 
IGM (|Hui fc Gnedinlll997l : Frheuns et al.ll2002h . In order to 
calculate it we assume that reionized gas is initially heated 
to 20, 000 K and subsequently evolves according to 



-A 



(3) 



where the number density of particles n changes due to adi- 
abatic expansion, and the right-hand side includes the vol- 
ume cooling rate A and heating rate Y. The cooling rate 
includes the contributions of single and dielectronic recom- 
binations, collisional excitations and ionizations, free-free 
emission, and Compton scattering off the CMB. Assuming 
a highly-ionized plasma in ionization equilibrium, the heat- 
ing rate for each atomic species equals its recombination 
rate times the average energy injected per photoionization. 
We assume an ionizing spectrum of stars with metallicity 
equal to 1/20 of the solar value, distrib uted according to 
the locally-measured IMF of lScalol l| 19981 ). This implies that 
helium is once-ionized simultaneously with hydrogen reion- 
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ization, and that later photoheating of recombining atoms 
results in an average energy injection of 3.8 eV per recombi- 
natio n to H I and 5.6 e V per r ecombination to He I [w here we 
used iLeitherer et all (|l999l )]. iHui fc Gnedlnl l| 19971 ) showed 
that after reionization, the gas in hydrodynamic simulations 
approximately follows a power-law temperature-density re- 
lation T oc p 7_1 at densities up to ~ 10 times the mean 
density. We estimate 7 by calculating the thermal evolution 
of gas at the (evolving) cosmic mean density, and the tem- 
perature of gas held at ten times the mean density. 

In general, star formation can be triggered when fresh 
gas is accreted from the IGM (including during mergers) or 
from merger-induced gravitational effects on gas that had 
previously accreted into galaxies. In this paper we assume 
that the SFR in galaxies is dominated by the fresh-gas mode, 
and comment below on the possibility of a significant mode 
of star formation in previous-accreted gas. We calculate the 
average accretion history of halos using the extended Press- 
Schechter model. Given that a halo of mass Mh formed at 
redshift Zh, the fraction of its mass that had already accreted 
into progenitor halos above Mmj n (z') at a given z' > Zh is 
l|Bond et al.lll99ll ; lLacev fe Cole] 1 19931 ) 



F acc (M h ,z h ,M min (z'),z') 



erfc 



S c (z') - 5 c (z h ) 



y/2 [S(M min (z')) - S(M h )] 



(4) 



where in general, S(M) is the variance of fluctuations in 
spheres of comoving radius containing mass M, and S c (z) 
is the critical density threshold required for a region to col- 
lapse by redshift z. We use the convention where the vari- 
ance is linearly-extrapolated to the present, while the critical 
threshold varies with redshift. If we write 8 c (z) = So(l + z), 
then So is approximately constant at high redshifts, and 
equals 1.28 for our assumed cosmological parameters. 

The SFH of a halo depends on the assumed reionization 
redshift z rc i of the surrounding region. For a halo containing 
a galaxy with stellar mass M* at Zh, the mean expected SFR 
at times earlier than reionization is simply equal to 



d 



SFR(z') = M t — F^ c (M h ,z h ,M mln {z'),z') 



(5) 



where t' is the age of the universe at redshift z' , and we as- 
sumed a star formation efficiency that does not change with 
time. We use this expression even though M m i n changes with 
time, since it only increases slowly with time and so we as- 
sume that gas that accretes into a halo near M m i n (z') at 
some z' is likely to continually merge into larger halos and 
thus be in a halo above M m i n (.z) at all z < z' . At reion- 
ization, however, the minimum mass for accreting new gas 
jumps suddenly to |Mj(z rc i), and a large fraction of the gas 
that had previously accreted remains for some time in halos 
less massive than this minimum infall mass. At a given red- 
shift 2' after reionization, the fraction of the halo gas that 
has accumulated since reionization equals the fraction that 
is in halos above jMj(z') and that had not yet accreted 
into galaxies by z TC i. For each possible progenitor mass at 
z' , we apply equation Q to find the fraction of the progen- 
itor's mass that still lay outside galactic halos at reioniza- 
tion. Note that this fraction is independent of the final halo 
mass Mh due to the ran dom walk statistic s of the extended 
Press-Schechter model l|Bond et al.lll99ll ). For progenitors 



of a given mass M at redshift z' , infall of new gas is due to 
both the increasing number of progenitors with time and the 
increase with time of the fraction of gas in the progenitors 
that did not already form stars before z la i. Thus, in this case 



SFR(z') = 



dM 



—F^M^z^M, z ) 



x {l-F acc (M, 2 ',Af min (z rci ),2 rci )} , (6) 

where M nl i n (z lc i) is calculated from eq. fTJ assuming either 
atomic or molecular cooling sets the minimum mass. 

Our approach in this paper to reionization is not to 
try to theoretically predict the absolute redshift at which 
it ends, since this depends strongly on the average escape 
fraction of photons from galaxies, a number that is highly 
uncertain. We instead focus on observational constraints on 
the reionization redshift as well as theoretical predictions of 
the spatial fluctuations in the reionization history. In par- 
ticular, we define a redshift shift Az which is the difference 
between the local redshift at which a given region is reion- 
ized and the redshift of complete global reionization. Regions 
in which galaxies form earlier than the cosmic average have 
a positive Az. The expected distribution of Az values can 
be predicted for a given reionization redshift, as long as the 
escape fraction from galaxies has roughly the same overall 
average value in different regions. We focus on these redshift 
fluctuations because they may be directly observable, as we 
argue below. 

Massive halos are rather rare at high redshift, but cur- 
rent telescopes can only detect the very brightest galaxies 
from those early times. At high redshift there are large 
fluctuations in the number densit y of galaxies in vari- 
ous regions l|Barkana fc Loebll20oi ). The extended Press- 
Schechter model predicts that galactic halos will fo rm ear- 
lier in regions that are overdense on large scales llKaisei 



1984l ; lBardeen et al.lll986l ; ICole fc K aiser 1989; Mo & White 



1996), since these regions already start out from an enhanced 



level of density, and small-scale modes need only supply the 
remaining perturbation necessary to reach the critical den- 
sity for collapse. Thus, an observed bright galaxy is likely 
to naturally flag a region with a large-scale overdensity. In 
such a region, galaxy formation and thus reionization will 
proceed faster than in the universe as a whole. In particu- 
lar, w ithin the extended Press-Schechter model l|Bond et al.l 
Il99ll ). if the cosmic fraction of gas accumulated in galaxies 
reaches a certain value at redshift z, then it will reach the 
same value at a different redshift z + Az in a box containing 
mass Mbox with a mean overde nsity Jb ox . At high reds hifts 
(z > 3), this shift in redshift is i|Barkana fc Loeb|[20o3 ) 



Az 



So 



(1 + z) x 



S(M b 



S{M min ) 



(7) 



To calculate the number density of halos above some 
mass in a given r egion, we begin w i th the cosmic mean halo 

~" dl999l) which fits numer- 



mass function of ISheth fc Tormen 



ical simulations mor e accurately 



Jenkins et ai1l200ll ) than 



the original model of lPress fc Schechteri (|l974l ). We then ad- 



just the halo distribution as a fu nction of density i n a re - 
gion based on the prescription of iBarkana fc Loebl l|2004) . 
in which the extended Press-Schechter model is used for the 
relative adjustment factor; we showed there that this hybrid 
model fits a wide range of simulations. If the expected num- 
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ber of halos above some mass in a given region is n, then 
the probability of finding at least one such halo in such a 
region is given by Poisson statistics: 1 ) = 1 — exp[— n]. 
Furthermore, since n is a function of the overdensity 5box of 
the region, the probability of the region having a given <5box 
can be calculated as the probability of having an overdensity 
<5box [which is given by a Gaussian with variance S(Mb ox )] 
times 1) for the given <5box, divided by P(^ 1) for the 
cosmic mean halo mass function. 

These probability distributions of Az also imply a 
large scatter in the reionization redshift in various regions. 
Once ionizing photons escape from the immediate vicin- 
ity of the sources, they encounter a mostly low-density 
IGM, where the dominant factor limiting reionization is the 
need to produce one ionizing photon per hydrogen at om 
IjMiralda-Escude. Haehnelt. fc Reej|200d : [Gnedinll200oh . If 
we assume a fixed efficiency for the production of ionizing 
photons in galactic halos, then the reionized fraction in any 
region is thus proportional to the total gas fraction in galax- 
ies within that region. This assumption is in fact fairly ac- 
curate even with recombinations included, as long as the 
number of times each hydrogen atom must be reionized is 
roughly the same in different regions. Indeed, the recombina- 
tion rate is not expected to vary greatly in different regions, 
since the fluctuation in the mean density of such regions is 
fairly small. For example, a 10 comoving Mpc box at red- 
shift 7 has density fluctuations of a ~ 0.2. Recombinations 
can be included in a future refinement of our model. 

Numerical simulations in small boxes i|Gnedinl |200Ch 
suggest that reionization occurs locally as a fast transition, 
as we have assumed for the progenitors of each galaxy, while 
the large-scale galaxy-number fluctuations mentioned above 
imply that glob ally reionization occurs at different times in 
different regions! Bar kana fc Loebl2004h . Rare over-dense re- 
gions can be assumed to reionize themselves as though they 
were isolated from the rest of the universe, since they are 
most likely to fully reionize before their lower-density sur- 
roundings. Thus, the Az calculated above yields an accurate 
estimate for the earlier reionization redshift of such regions, 
as long as Az is positive and significantly large. We assume 
this is the case for regions that are at least a +1 — a fluctu- 
ation, which means that they form more galaxies than 84% 
of other such regions in the universe. For voids, on the other 
hand, while Az of reionization can come out negative in our 
model, the real Az in the universe cannot be negative for 
any region; as soon as the universe contains enough ionizing 
photons to fully reionize, all regions are indeed reionized, 
including the underdense regions which are reionized by ex- 
cess photons from the denser surrounding regions. Thus, we 
adjust the probability distribution of Az for reionization, for 
all Az smaller than that corresponding to a +1 — a fluctu- 
ation, by a constant factor so that the total probability of 
having Az > is unity. 

These arguments suggest that the SFHs of high-redshift 
galaxies can be used to study the spatial correlation func- 
tion of reionization. In a sufficiently small region, all galax- 
ies will show in their SFH evidence for approximately the 
same reionization redshift which characterizes that partic- 
ular region. Galaxies separated by a large distance will in 
general show a significant variation in reionization redshifts, 
where this variation will be strongly correlated with the dif- 
ferent environments of the galaxies. In the quantitative re- 



sults below, we adopt a distance of 10 comoving Mpc as the 
fiducial size of a region with a uniform reionization history, 
based on analytical m odels of the distribution of bubble sizes 



during reionization ( 


Furlanetto. Zaldarriaea, & Hernauistl 


l2004lWvithe & Loeb! 


2004V 



3 RESULTS 

The reionization signature we predict is a generic feature, 
and we illustrate it in this section for a wide range of possi- 
ble galaxies. For example, Figure[T]shows the mean expected 
SFHs of redshift 6.5 galaxies. As detailed in the previous 
section, we assume that stars form in halos above a mini- 
mum cooling mass before reionization, and in halos above a 
minimum infall mass after reionization. We consider a pre- 
reionization minimum mass set by requiring atomic cool- 
ing or molecular cooling. We find that a massive z — 6.5 
galaxy should generically show a double-peaked SFH, with 
one early peak due to star formation in low-mass progenitors 
before reionization, and a much later peak due to gas infall 
into more massive progenitors. The same qualitative fea- 
tures are expected even when we vary the observed redshift 
(sec Figure [5J) or the reionization redshift (see the bottom 
panel of Figure [3} . For a given observed galaxy, if reioniza- 
tion occurred earlier, then its signature on the SFH is more 
difficult to observe, for two reasons. First, a smaller fraction 
of the observed stars is expected to have formed prior to 
reionization (see the top panel of Figure [3}; and second, the 
early peak in this case lies more than a few hundred million 
years in the past, allowing only low-mass stars to survive 
and making it difficult to pin down the age of this episode 
of star formation. 

Current instruments are limited to detecting the bright- 
est high-redshift galaxies, which should typically reside in 
the most massive halos. Such halos are far more likely to 
form in overdense regions, where galaxy formation generally, 
and reionization in particular, should occur earlier than in 
the universe as a whole. Figure U shows that substantial val- 
ues of Az up to ~ 7 are possible, depending on the halo mass 
and the redshift at which it is observed. Note that the value 
of Az is essentially independent of the reionization redshift 
and of the precise value of the minimum mass, since for 
both atomic and molecular cooling, S(M nl i n ) S> •S'(Mbox) 
in eq. (0). Increasingly rare halos should show larger val- 
ues of Az, a trend that should be seen observationally once 
large-area surveys determine the abundance of rare, mas- 
sive halos. Figure [5] shows the mean abu ndance of halos at 
several redshifts l|Sheth fc Tormenlll999h . 



4 APPLICATION TO PARTICULAR 
GALAXIES AND DISCUSSION 

iMobasher et al.l (|2005l ) have recently used deep multi- 
waveband photometry to place the galaxy HUDF-JD2 po- 
tentially at a redshift z ~ 6.5. They estimated a stellar mass 
of 6 x 10 Mq , for a Salpeter stellar IMF which places most 
of the mass in stars below a solar mass. Fixing the con- 
tribution of the 1-3Mq stars which dominate the observed 
in frared emiss ion, we convert to the locally-measured IMF 
of lScalol (|l998h and find a total stellar mass of 3 x 10 n M Q . 
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Figure 1. SFHs of redshift 6.5 galaxies vs. cosmic age or red- 
shift. We show the fractional SFR summed over all progenitors 
that form stars. Assuming 2 rc ; = 10, we consider atomic cool- 
ing for the progenitors of = 10 13 Mq halos (solid curve) or 
iWb = 10 11 A/q halos (dotted curve), or molecular cooling for 
the progenitors of = 10 13 Mq halos (long-dashed curve) or 
Mjj = 10 11 Mq halos (short-dashed curve). 
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Figure 2. Same as Figurc[T] but for redshift 8 galaxies (top panel) 
and for redshift 5.8 galaxies (bottom panel). Note that in each 
panel the x-axis shows cosmic age (bottom) and redshift (top). 



Figure 3. Bottom panel: Same as Figures [pj and [2] (with the 
x-axis shown as both time and redshift), but for the case of a 
10 13 A/q halo and atomic cooling, assuming no reionization (solid 
curve) or 2 re ; = 9, 10, or 13 (dashed curves). We show the mean 
expected SFR in absolute units, assuming a 20% star formation 
efficiency. Top panel: Fraction of stars that formed prior to reion- 
ization vs. 2 re ;, for a halo observed at z = 6.5. We consider atomic 
cooling for halos M h = 10 13 A/ Q (solid curve) or M h = 10 11 M Q 
(dotted curve), and molecular cooling also for = 1O 13 M0 
(long-dashed curve) or Mj, = IO^Mq (short-dashed curve). 



This implies a minimum host halo of total mass 2 x 10 1 M@, 
for a maximum possible (100%) efficiency of formation of 
stars from baryons. Efficiencies near ~ 10% are more typ- 
ically assumed, so we illustrate results f or an intermediate 
value of 20% which yields a 10 13 A'/ Q halo. lEvles et al.l (|2005l ) 
have also analyzed two high-redshift galaxies with spectro- 
scopic redshifts z = 5.78 and z = 5.83, with stellar masses 
> 10 10 A/ Q [for the lScalol <|l998D IMF], which for a 20% effi- 
ciency implies host halos of mass 3 x 10 n MQ for each galaxy. 

Galaxies similar to those observed should generally 
show a clear double-peaked SFH. For a 1O 13 M0 galactic 
halo at z — 6.5, a majority of the stars are expected to 
have formed prior to reionization if reionization occurred 
at z re j < 9.9 or z rc i < 11.3, assuming atomic or molecu- 
lar cooling, respectively (see Figure These redshifts cor- 
respond to look-back times of 363 and 443 million years, 
respectively, with respect to the cosmic age of 847 million 
years at z = 6.5. A 3 x lO n M halo observed at z = 5.8 
should have its SFH dominated by pre-reionization stars if 
Zrci < 8.9 (atomic cooling) or z Ic i < 10.4 (molecular cool- 
ing). Indeed, signatures of e arly star formation were found 
in the galaxies cited above. lEvles et al l {2005) found that 
most of the stars in the two z = 5.8 galax ies had formed 
at z > 7.5. Similarly. iMobasher et a l. (2005) found that the 
z — 6.5 galaxy had formed most of its stars at z > 9, and 
could h ave reionized its sur rounding region even earlier, at 
z > 10 (|Panagia et al.ll2005l ). 

If these signatures do indeed correspond to local reion- 



6 



R. Barkana and A. Loeb 




Figure 4. Bottom panel: Most probable redshift shift vs. halo 
mass, for a halo observed at z = 5.8 (long-dashed curves), 6.5 
(solid curves), or 8 (short-dashed curves). In each case, we show 
Az for reionization of a 10 Mpc box containing the halo, relative 
to complete reionization of the universe. Top panel: Normalized 
probability distribution of Az. We consider halos = 10 13 Mq 
(solid curves) or M h = 10 12 Mq (dotted curves), considering 
(from left to right) a halo observed at z = 6.5 or at z = 8. For the 
Mjj = 10 11 Mq case, we also consider z = 6.5 (short-dashed) or 
z = 8 (long-dashed). Note that to the left of the dotted vertical 
line (which indicates a +1 — a fluctuation), the total probability 
in each case is fixed but its precise distribution is uncertain (see 
§ 2). 

ization of the regions surrounding these galaxies, then they 
could have important implications for cosmic reionization. 
However, in analyzing these implications we must account 
for the fact that the galaxies were most likely found in over- 
dense regions that reionized earlier than did the rest of the 
universe. Indeed, the 10 comoving Mpc region surround- 
ing the massive z — 6.5 galaxy most likely reionized at a 
Az = 4.7 earlier than the rest of the universe, i.e., at z > 11. 
The observed z — 5.8 galaxies correspond to halos that are 
not nearly as rare, but still are most likely in regions that 
reionized early by a Az — 1.2, i.e., at z > 7.5. Thus, the sig- 
natures of reionization observed in these galaxies only imply 
an early reionization of the surrounding regions, but do not 
push cosmic reionization earlier than the already established 
z > 6.5 (see § 1). 

If these galaxies indeed lie in atypical, overdense re- 
gions, then consistency implies that such massive galaxies 
must be extremely rare. The two z = 5.8 galaxies were 
detected i n an effective volume of ~ 1.7 x 10 comov- 
ing Mpc 3 (jStanwav, Bunker, fc McMahonll2003l ). while the 
z = 6.5 galaxy was found in the Hubble Ultra Deep Field, 
which (adopting an effective redshift range of 6.2-6.8 for the 
search) corresponds to ~ 9 x 10 3 comoving Mpc 3 . Although 
the number of galaxies is small and large fluctuations are 
expected, we can estimate the resulting abundances of these 
galaxies (see Figure [5}. The z = 5.8 galaxies have a lower 



Figure 5. Cosmic mean number density of halos versus halo 
mass. We show the comoving number density at z = 5.8 (dotted 
curve), 6.5 (solid curve), 8 (short-dashed curve), 10 (long-dashed 
curve), and 15 (dot-dashed curve). Also show n are the observed 
abundances of massive z = 5.8 galaxies using Eylcs et al. (2005) 
and assuming a 20% star formation efficiency (x), and the range 
(two connected circles) corresponding t o an efficiency of 20% - 
100% for the massive z = 6.5 galaxy from lMobasher et al.l J2005h ■ 



abundance than expected, by a factor of ~ 20. However, the 
observed number of two galaxies represents a lower limit, 
since there are other galaxies in the field with undetermined 
redshifts or stellar masses. Furthermore, if the star forma- 
tion efficiency in these galaxies is only ~ 5% rather than 
20%, then their abundance is similar to that expected. How- 
ever, the z = 6.5 galaxy is extremely massive and was found 
in a very small field. For the 20% efficiency that we have 
adopted, the a posteriori probability of finding this galaxy 
is ~ 6x 10 -6 , making it rather unlikely that the claimed high 
redshift is accurate. If this galaxy represents the result of an 
unusual star formation episode which was characterized by 
a near-perfect efficiency, then finding such a 2 x 1O 12 M0 
halo represents a ~ 10~ 3 likelihood. Further observations of 
larger fields are required to determine whether this detec- 
tion may have been a stroke of luck. Note, however, that 
even a small change in the halo mass would change the de- 
tection probability significantly, because of the exponential 
cutoff of the halo mass function. Thus, the probability is also 
sensitive to any uncertainties in the estimated stellar mass. 

We caution that the SFHs we predict in this paper may 
not be entirely accurate. First, even in a small region from 
which a single observed galaxy arises, the suppression of star 
formation builds up somewhat gradually as the progenitors 
reionize their surroundings. In addition, if significant star 
formation is triggered during mergers in previously-accreted 
gas then this will somewhat weaken the signature of reion- 
ization, since infall of fresh gas is immediately suppressed 
by the increased pressure while mergers continue among 
progenitors that had accreted gas prior to reionization. We 
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expect, however, that infall will be dominant at high red- 
shift, since around reionization most of the gas of an even- 
tual large galaxy lies in low-mass progenitor halos (typically 
with mass of order lf/M©) in which even a few supernovae 
can easily disperse much of the gas. Thus, although more 
detailed models of feedback should be explored, we expect 
the initial burst of star formation in fresh infalling gas to 
be the dominant mode of star formation at high redshift. 
It is also important to note that we have calculated in this 
paper mean expected SFHs, while there should be signifi- 
cant variation in the merger histories of individual galaxies 
(e.g., see the top panel of Fig. [4]). Detailed models and nu- 
merical simulations can be used to more accurately quantify 
our predictions. However, regardless of the fine details, we 
emphasize that the signature we predict is a generic, quali- 
tative feature that should be clearly observable in a majority 
of high-redshift galaxies. 

In conclusion, observations of high-redshift galaxies 
should always be fitted with SFHs consisting of at least two 
separate episodes of star formation. If the signature of reion- 
ization is seen in the SFHs of many individual galaxies, then 
the variation of the reionization redshift in different regions 
and its spatial correlation function can be studied quanti- 
tatively. In particular, a consistent picture should emerge 
where overdense regions both reionize earlier and contain a 
higher abundance of massive galaxies than the rest of the 
universe. 
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